Speckle-visibility spectroscopy: A tool to study time-varying dynamics 
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We describe a multispeckle dynamic light scattering technique capable of resolving the motion of 
scattering sites in cases that this motion changes systematically with time. The method is based on 
the visibility of the speckle pattern formed by the scattered light as detected by a single exposure 
of a digital camera. Whereas previous multispeckle methods rely on correlations between images, 
here the connection with scattering site dynamics is made more simply in terms of the variance of 
intensity among the pixels of the camera for the specified exposure duration. The essence is that 
the speckle pattern is more visible, i.e. the variance of detected intensity levels is greater, when 
the dynamics of the scattering site motion is slow compared to the exposure time of the camera. 
The theory for analyzing the moments of the spatial intensity distribution in terms of the electric 
field autocorrelation is presented. It is demonstrated for two well-understood samples, a colloidal 
suspension of Brownian particles and a coarsening foam, where the dynamics can be treated as 
stationary. However, the method is particularly appropriate for samples in which the dynamics 
vary with time, either slowly or rapidly, limited only by the exposure time fidelity of the camera. 
Potential applications range from soft-glassy materials, to granular avalanches, to flowmetry of living 
tissue. 
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I. INTRODUCTION 

Dynamic light scattering (DLS) is a powerful tool for 
probing motion within samples of physical, chemical, bi- 
ological, and medical interest Q, 0, fi S II B ■ The 
physical basis is that the frequency spectrum of the scat- 
tered light is Doppler broadened according to the veloci- 
ties of all the scattering sites. The shape of the spectrum 
reveals the nature of the motion, for example whether 
it is ballistic or diffusive; the characteristic width of the 
spectrum reveals the rate of the motion, for example the 
root-mean-squared speed or the diffusion coefficient. If 
the sample is nearly transparent, so that incident photons 
scatter at most once, then the spectrum can be resolved 
vs scattering angle in order to probe collective motion 
at different length scales. This is the single-scattering 
regime. By contrast if the sample is opaque, so that in- 
cident photons scatter off many sites before exiting the 
sample, then any wavevector-dependence is lost. The art 
of DLS in this re gim e is known as diffusing-wave spec- 
troscopy HHHulu). 

The most straightforward approach to DLS is to mea- 
sure the frequency spectrum directly, for example using 
a Fabry-Perot interferometer with a very narrow band 
pass. However, it is also common practice to deduce 
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the spectrum by an interference technique, in which the 
scattered light is collected over an area comparable to 
one speckle spot (spatial-coherence length). The motion 
of the scattering sites causes corresponding changes in 
the speckle pattern, and hence large fluctuations in the 
detected intensity. These fluctuations are quantified by 
the temporal intensity autocorrelation function, which 
is simply related to the frequency spectrum under cer- 
tain conditions (below). This is known as intensity- or 
photon-correlation spectroscopy (PCS). One advantage 
of PCS is that digital correlators are commercially avail- 
able that can compute the intensity autocorrelation over 
many decades in delay time, for example 10 ns to 100 s. 
One disadvantage of this approach is that the tempo- 
ral fidelity is limited by the necessity of sampling over 
many correlation times to build up statistical weight. 
This makes them a poor choice for studying systems with 
dynamics changing on time scales of seconds or faster. In- 
terferometers are useful for large frequency shifts, but do 
not sport an equally impressive dynamic range. Given 
the breadth of applications of DLS, it is perhaps not 
surprising that the essential equivalence of information 
available from interferometric and correlation-based ap- 
proaches to DLS is not universally recognized 01 • 

In order to ensure simple connection between the in- 
tensity autocorrelation and the frequency spectrum of 
the scattered light, several conditions must be met: (a) 
There must be many, uncorrelated scattering sites or re- 
gions; (b) the extent of the motion must be sufficiently 
great as to fully randomize the speckle pattern; and (c) 
the scattering site dynamics must not vary over the time 
scale of the measurement. The first criterion holds if the 
sample and scattering volume are sufficiently large; this 
does not represent a fundamental restriction. The second 
criterion holds if the sample is fluid or if the scattering 
sites are bound only loosely to a fixed average location. 
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The third criterion holds if the sample is in thermal equi- 
librium, or if the sample is in a stationary state in which 
both the external energy input and the microscopic dy- 
namical response do not fluctuate. Thus, the conditions 
(a-c) for conventional single-detector PCS to apply are 
not overly restrictive. It is possible to verify whether not 
these conditions hold through measurement of higher- 
order temporal intensity correlations [13j, which can be 
processed from the raw intensity vs time data stream 
simultaneously with the second-order intensity autocor- 
relation. 

There are many systems where some of the above con- 
ditions do not hold and conventional single-detector PCS 
does not apply. The kinetics of phase separation, gela- 
tion, and aggregation are examples of long-standing in- 
terest, in which the dynamics progressively change with 
time [l4|. These processes can be treated as stationary 
only if the evolution is slow compared to the time scale 
over which the intensity autocorrelation decays. The 
broad class of soft-glassy materials comprise another ex- 
ample where the dynamics change with time Fur- 
thermore, just as for the gelation problem, the scattering 
sites can become more tightly bound with age, so that 
after a certain point the speckle pattern no longer fully 
randomizes. And lastly, dynamics in granular materials 
usually cannot be studied with traditional PCS, for ex- 
ample because the input of energy is vibratory or because 
the response is intermittent avalanche-like flows 

Multispeckle dynamic light scattering techniques have 
been introduced as a useful remedy in such situations 
where traditional single-detector PCS methods do not 
apply [13 list - The approach is to compute the tem- 
poral autocorrelation function for each pixel of a digital 
camera, and then to average together the results. Since 
there are many pixels, and hence many speckles, it is 
no longer a requirement that motion within the sample 
cause the speckle pattern to fully randomize. And since 
the large number of pixels can significantly reduce the 
time needed to acquire good signal-to-noise, it is easier 
to study evolving dynamics. However it still remains a 
challenge to implement multispeckle DLS. A prohibitive 
difficulty is that commercial multispeckle devices do not 
exist. A limiting difficulty is either that vast quantities 
of data must be stored for post-processing or that real- 
time processing must be made sufficiently fast. Further 
difficulties arise from the use of charge coupled devices 
as light sensing elements. Hardware and software ad- 
vances continue to be reported in the technical literature 

[S El IH El El El M H|. 

In this paper we supply full details and demonstra- 
tion of a multispeckle dynamic light scattering technique 
we dub speckle-variance spectroscopy E3 or speckle- 
visibility spectroscopy (SVS) Ell- Our approach is to 
characterize motion within a sample in terms of the vis- 
ibility of the speckle pattern formed with scattered light 
for a single exposure of a digital CCD or CMOS cam- 
era. We begin by introducing appropriate notation and 
the experimental apparatus in the context of the more 



usual multispeckle DLS. Then we describe the theoreti- 
cal underpinnings of SVS, and give examples for common 
types of scattering site dynamics. Our theory contra- 
dicts a widely-cited prediction obtained in the context 
of laser-speckle photography E3- Next, crucially, we 
demonstrate the validity of our theory by experiments 
on well-known samples. Finally we discuss experimental 
considerations for successful implementation of dynamic 
light scattering with a digital camera. 



II. PHOTON-CORRELATION 
SPECTROSCOPIES 

We begin with prerequisite theoretical and experimen- 
tal background materials necessary for the next sections 
on SVS. 



A. Theory 

In all DLS experiments, light from a coherent source 
enters a sample. Some portion scatters, with individ- 
ual photons experiencing different trajectories, and some 
fraction of the scattered light reaches a photodetector. 
Ignoring constant factors, the detector reports a signal 
proportional to the light intensity, I(t) — E(t)E*(t), 
where the electric field E(t) is a superposition of many 
fields representing many photon trajectories. The ac- 
quired intensity can be an analog signal, or it can be a 
bit-stream with each pulse representing a different de- 
tected photon. Ultimately, the quantity of interest is 
either the power spectrum |i?(ti;)| 2 or its Fourier trans- 
form: the temporal electric field autocorrelation function. 
We denote the absolute normalized temporal electric field 
autocorrelation as 



9i{t) = \(E(t)E*(t + T))\/(E(t)E*(t)), 



(1) 



where r is the delay time. In traditional PCS the aver- 
age (• • •) is taken over a range of times t stai t < t < t st0 p- 
By definition, gi(r) decays from one to zero as r ranges 
from zero to infinity. The characteristic time scale for 
the decay is the reciprocal of the characteristic width of 
the power spectrum. If the power spectrum is symmet- 
ric and centered around u) a , then the normalized (but 
not absolute) electric field autocorrelation function is 
g\(r)e lM ° T . For example, a Lorentzian power spectrum 
|-E(w)| 2 ex l/{Dq 2 + (uj— ll>o) 2 ] and an exponential field au- 
tocorrelation function <7i(r) = exp(— Dq 2 r) correspond 
to light of incident frequency uj q scattered by wavevec- 
tor q from diffusing particles; the diffusion coefficient D 
could be extracted from measurement of either the power 
spectrum or the field autocorrelation. 

In single-detector PCS the electric field autocorrelation 
function is deduced from measurement of the normalized 
intensity autocorrelation function, 



g 2 (T) = (I(t)I(t + T))/(iy< 



(2) 
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This is straightforward only if the three conditions dis- 
cussed in the Introduction are all met. If (a) the elec- 
tric field is the superposition of many independent scat- 
tered fields and if (b) the field autocorrelation decays 
to zero over a time scale much shorter than the dura- 
tion of the measurement, then the Central Limit The- 
orem implies that E(t) is a Gaussian-distributed com- 
plex variable with zero mean. Intuitively, the total field 
E(t) = at some instant of time, t, may be 

evaluated graphically by phasor addition. If there are 
enough independent scattering regions, then each term 
in the sum constitutes one step in a random walk in 
the complex plane. Many such random walks will be 
sampled, and hence the distribution of values of E{t) 
over the course of the measurement will be Gaussian, 
if <7i(t) fully decays to zero over a shorter time scale 
than the measurement duration. If the field distribu- 
tion is Gaussian, then temporal correlations of the form 
(E(t)E*{t)E{t + T 1 )E*{t + T 1 )E{t + T 2 )E*(t + T 2 ) ■ • •) can 
be expressed as a sum of products of field autocorrela- 
tions. For example, the normalized intensity autocorre- 
lation is a four-order field correlation that reduces to 



32 (r) = l + (r)] 2 , (3) 



where /3 < 1 is a number determined by the ratio of detec- 
tor size to speckle spot size. This is widely known as the 
Siegert relation. A detailed derivation of the Siegert rela- 
tion, the value of /3, and analogous results for third- and 
fourth-order temporal intensity correlations, are given in 
Ref. To briefly summarize, the method of PCS is to 
measure gi (r) and to extract g\ (r) using Eq. @ . Subse- 
quent connection is then to be made between g\ (r) and 
scattering site dynamics, depending on details of the il- 
lumination and detection geometry and on the optical 
properties of the sample. 

In multispeckle PCS methods the intensity autocorre- 
lation is measured at each pixel of a digital camera and 
the results are averaged together [H EH E E3, EJ |H 
I231 0, |2|| . By virtue of the large number of pixels, the 
combined statistics of all the detected fields is now Gaus- 
sian even if the field autocorrelation never decays to zero. 
In effect, statistics are sampled by an ensemble average 
over many speckles rather than by a time-average for a 
single speckle. Thus the Siegert relation, Eq. ©, may be 
invoked even more generally for multispeckle PCS than 
for single-detector PCS. 

As an aside, the violation of the Siegert relation in 
single-detector PCS due to non-randomization of the de- 
tected electric field is sometimes said to be due to non- 
ergodicity of the sample. This is a misnomer and can 
lead to confusion. The ergodicity of dynamics within the 
sample, and the ergodicity of the field statistics for the 
detected light, are distinct issues that may or not be re- 
lated. 



B. Experiment 

We now apply the above multispeckle PCS technique 
to a suspension of diffusing Brownian particles. This 
serves as a starting point from which to demonstrate 
SVS, since all measurement and sample hardware carry 
over without change. 

The sample consists of 653 nm diameter polystyrene 
spheres [Duke Scientific] suspended in water at a vol- 
ume fraction of ten percent. It is poured into a glass 
beaker, diameter 6 cm, to a depth of 2.4 cm, then sealed. 
Light from a Coherent Verdi V5 NdYV04 Laser, wave- 
length A = 532 nm, is expanded and directed almost 
normal to the bottom of the sample beaker with a Gaus- 
sian spot size of a = 1.25 cm. The outpower of the 
laser is held fixed, and the illumination intensity is re- 
duced as needed by use of neutral density filters. See the 
schematic diagram in Fig. ^ According to Mie scatter- 
ing theory for dilute independent spheres [3(]|, the scat- 
tering length specifying the exponential attenuation of 
a beam is l s = 24 /im, and the average cosine of the 
scattering angle is g — 0.90. Therefore, about ten scat- 
tering events are required to randomize the photon prop- 
agation direction, and the transport mean free path is 
I* = 's/(l — g) ~ 240 /im. Thus our sample has an 
opaque white appearance, and we operate in the multiple 
scattering regime known as diffusing- wave spectroscopy. 

In order to perform multispeckle DWS, a portion of 
the backscattered light leaving the bottom of the sample 
is reflected by mirror into a Basler-160 digital line scan 
CCD camera. This device has 1024 pixels, each 10 /im 
x 10 /im and 8 bits deep, and can capture images at a 
maximum rate of 58 kHz. Except for the mirror and a 
532 nm line filter, there are no other optics. The sample- 
to-camera distance is adjusted to about d = 30 cm. This 
gives a speckle size of s k dX/a = 13 /im, and a ratio 
of pixel to speckle areas of A p - lxe i/ ^4 spec ki e ~ 0.6. The 
camera is interfaced to a PC equipped with a National 
Instruments PCI-1422 card, and is programmed using 
LabVIEW. 

As a benchmark reference to compare with our SVS 
technique, our operating procedure is to record the inten- 
sity levels in all 1024 pixels for a total of 2 s in increments 
of 20 /is; the entire data set thus consists of 102.4 million 
8-bit values. The laser intensity is adjusted so that the 
average gray-scale value is 40. When the laser is blocked, 
the signal drops to a "dark count" grayscale value of 3.5. 
The first step in the analysis is to subtract the dark count 
and divide by the average remaining signal, thus giv- 
ing the normalized intensity time trace I(t) / (I) for each 
pixel. A portion of such a trace for one pixel is displayed 
in Fig. |21 As the colloidal particles diffuse, the intensity 
level at a pixel indeed fluctuate strongly; here, it is seen 
to vary between about 0.4 and 4 times the average. 

The observed intensity fluctuations in Fig. |2] display 
structures lasting over a range of time scales. As is done 
in multispeckle PCS, this behavior may be quantified by 
the normalized intensity autocorrelation, gi (t) defined 
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FIG. 1: Experimental apparatus for performing both multi- 
speckle dynamic light scattering and speckle-visibility spec- 
troscopy on an opaque colloidal suspension. The orientation 
of the array of pixels is horizontal, perpendicular to the plane 
of the drawing. 
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FIG. 2: An example of intensity vs time in one pixel (left) and 
the normalized intensity autocorrelation function | pi (t) | 2 = 
[(7(0)7(r)}/(/) 2 - (right), for light backscattered from 
an opaque colloidal suspension. The intensity autocorrela- 
tion is directly computed from the full time-trace, which con- 
sists of 10 points in 20 fis time increments for each of the 
1024 pixels in the linescan CCD camera. The normalization 
factor /3 = 0.34 was obtained by extrapolating the unnormal- 
ized intensity autocorrelation to zero delay time. The sus- 
pension consists of 653 nm diameter polystyrene spheres in 
water, at a volume fraction of 10 percent. The theory of 
diffusing-wave spectroscopy predicts <?i(t) ~ exp[— 7^/ 6t/t ] 
where t = 1/Dk 2 =6.1 ms. The fit to this form (dashed 
curve) gives 7 = 1.3. 



by Eq. which we compute directly for each pixel and 
then average together. According to the Siegert relation, 
Eq. ©, the zero-time intercept is (72(0) = 1 + (3. Extrap- 
olating <?2(t) data to r = gives f3 = 0.34, which is con- 
sistent with the ratio of pixel to speckle areas. Invoking 
the Siegert relation, we deduce the normalized field auto- 
correlation |<7i(t)|, and plot its square in Fig. [2 Finding 
the value of j3 is often referred to as the issue of nor- 
malization. Note that the square of <?i(t) is simply the 
intensity autocorrelation displayed dimensionlessly on a 



scale ranging from 1 to 0. Thus, the time scales of struc- 
ture in the intensity time trace can be compared directly 
with features in the decay of |<7i(t)| 2 . Indeed the decay is 
very fast initially, reflecting the fast fluctuations in I(t). 
The later-time decay is slower, reflecting the longer-lived 
fluctuations evident in Fig. [21 as drift in a local average 
of the intensity. 

The above measurement of gi(r) may be compared 
with the predictions of DWS. In the backscattering geom- 
etry, with equivalent plane-wave in / plane-wave out illu- 
mination and detection, the theory of DWS predicts 
gi(r) rj exp(— 7\/6t/t ) where typically 1.5 < 7 < 2.5 
and where r = l/(Dk 2 ) is the characteristic time for a 
particle to diffuse a distance 1/k — A/(27rn) where n is 
the index of refraction. For our sample, the predicted de- 
cay time is r =6.1 ms. The stretched-exponential form 
of gi (t) reflects the broad length distribution of possible 
photon paths that contribute to the signal. It also reflects 
a subtle breakdown of diffusion and continuum approx- 
imations for short path lengths [3ll l32|| . The value of 
7, but not the stretched-exponential form, is particularly 
sensitive to the treatment of short paths and can thus 
be affected by the polarization states, boundary reflec- 
tivities, and propagation directions for the incoming and 
outgoing photons [331,03 113 Taking the value 7 = 1.28, 
somewhat lower than expectation, we obtain an excellent 
fit to the data as shown by the dashed curve in Fig. [2] 
Thus we pronounce our sample, apparatus, dataset, and 
multispeckle analysis methods as sound. For demonstra- 
tion of SVS in a later section, the value of 7 will not be 
important; we only need a sample with known gi (r). 



III. SPECKLE- VISIBILITY SPECTROSCOPY 

In this section we develop the theory of speckle- 
visibility spectroscopy (SVS). The underlying principle 
of SVS is illustrated in Fig. [31 which displays intensity vs 
pixel number for four different exposure times, T, of the 
camera. The shortest exposure in Fig. [21 T — 2 x 10~ 5 s, 
is shorter than the decay time of gi (r); therefore, the 
speckle appears static and large intensity differences are 
registered from pixel to pixel. For longer and longer expo- 
sures, the visibility of this speckle pattern progressively 
fades. This is because the intensity at each individual 
pixel fluctuates during the exposure and is averaged over 
the exposure time. In the limit of a very long exposure 
time in comparison to the decay time, each pixel ap- 
proaches the same mean intensity value and there is no 
variation among the pixels. Indeed, the longest exposure 
in Fig. [21 T = 2 x 10~ 2 s, is longer than the decay time of 
<7i(t); here, many speckles are sampled at each pixel over 
the duration of the exposure, and each pixel registers a 
value close to the average. The very essence of SVS is, 
thus, to quantify the visibility of the speckle pattern in 
terms of moments of the distribution of intensities regis- 
tered by all the pixels for a given exposure duration, and 
to relate this to the absolute normalized electric field au- 
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FIG. 3: Intensity vs pixel number, ie the profile of the speckle 
pattern in the plane of the CCD camera, for the same colloidal 
suspension and optical configuration as in Figure 1. The expo- 
sure durations T differ by successive factors of ten, as labelled. 
Since the speckles change with time, as shown in Fig. [5] their 
visibility is smaller for longer exposures. This is the essence 
of SVS. 



tocorrelation function gi(r). A subsequent connection 
with scattering site dynamics can then be made as per 
usual DLS practice in either single- or multiple scattering 
limits. 

Before carrying out the theoretical aspects of this pro- 
gram, we note that our method is not without precedent. 
Perhaps the first is a calculation |3Jj and experimental 
verification |37| of the distribution for the photocurrent 
as measured by one detector as a function of integra- 
tion time. Another precedent is "laser-speckle photog- 
raphy" [2^|, in which the blurring of speckle in a laser- 
illuminated scene is taken as a signature of motion 0, Q • 
The latter is now being applied to cerebral blood flow, 
in particular [3^, |3^, 0, |4l|. One aspect of our con- 
tribution here is to simplify and generalize the work of 
Refs. |3^, and to correct a mistake in the widely-cited 
work of Ref. El. 



A. Variance 

The variance of intensity across the pixels is a sim- 
ple way to quantify the visibility of the speckle pattern 
formed at the imaging array. For a given exposure, each 
pixel reports a signal that is proportional to the total 
number of photons it receives. Thus the signal at pixel i 
is proportional to the time-average of the intensity trace 
/,•(*): 



SiT — 



Ii(t')dt'/T, 



(4) 



where t = defines the beginning of the exposure and T 
is the duration of the exposure. The data returned by the 
camera, for a single exposure, consists of the set {S^t} 
where the index i ranges from 1 to the total number N of 



pixels. All quantities of interest are to be computed from 
the N members of this set. For example the n -moment 
of the distribution of pixel signals is 



(5) 



where the subscript T is a reminder that the result de- 
pends on the exposure duration. Note that these mo- 
ments represent an ensemble average over pixels for a 
fixed time interval. 

To compute the variance we focus on the first two mo- 
ments of the signal distribution. The first moment is 
simply the average intensity, (I) — J2iLi S%,t/N, which 
is independent of the exposure duration. The second mo- 
ment is the average over pixels of the quantity 



(s hT ) 2 = r r i i (t , )i i (t")dt'dt"/T 2 

Jo Jo 



(6) 



Since this is an ensemble average, the Siegert relation 
Eq. © may be invoked: (ii (*')£(*")) = {I) 2 {l+f3[g l {t l - 
t")] 2 }, giving an intermediate result for the second mo- 
ment as 



(If 



{l+(3[ gi (t'-t")} 2 }dt'dt"/T 2 . (7) 



The first term in the integral is one; the second term can 
be reduced to a single integral by recognizing that gi(t) 
is usually an even function. We now define a normalized 
variance, and finish the calculation: 



V 2 (T) 



\ W 2 )tI(i? 1] 

T i-T 



I I [ 9 i(t'-t")] 2 dt'dt"/T 2 , 
Jo Jo 



2(l-t/T)[ 9l (t)] 2 dt/T. 



(8) 



This is the fundamental equation of SVS. The top line is 
a definition; it quantifies speckle visibility on a scale of 
— 1 in terms of the first two moments of the distribution 
of pixel signal data, {S^y}, returned for a given exposure 
of duration T. The middle line is an intermediate step 
that holds even if gi(t) is not even. The bottom line is 
where contact usually is to be made between measure- 
ment and the underlying normalized electric field auto- 
correlation. Evidently the variance is a weighted-average 
of [gi (t)] 2 over the exposure interval < t < T, with 
heavier weighting for shorter t. This weighting reflects 
the distribution of possible time differences within an ex- 
posure. 



B. Higher-Order Moments 

The distribution of pixel signals is typically skewed to- 
ward higher values, as seen for example in Fig. |3J there- 
fore, it is not Gaussian and cannot be fully specified by 



6 



just the value of the variance. Hence we repeat the cal- 
culation leading to the fundamental equation of SVS, 
Eq. JHJl, but now for higher-order moments. The results 
can be useful for diagnosing problems with the experi- 
mental apparatus, for deducing the normalization con- 
stant (3, and for better testing trial forms of gi(t) for 
unknown samples. 

We define reduced moments of the pixel signal distri- 
bution as 

v n (T) = (I n )r/(I) n 1. (9) 



The reduced second moment is recognized as i)2(T) — 
(3V% (T) . The third and fourth reduced moments, by con- 
trast, contain terms that are not proportional to (3\ there- 
fore, their dependence on the number of speckles per pixel 
cannot be normalized away by a simple division. 

Note that for very short exposure times, the reduced 
moments approach 

v n (0) = [1 + ft [1 + 2/3] • • • [1 + (n - 1)13] - 1, (12) 

which is a well-known result describing the moments of 
static speckle patterns F° r l° n S exposure times, 

the reduced moments vanish as 

/>oo 

v 2 (T) -> 2(3 J [ 9l (t)} 2 dt/T, 
v 3 (T) -> 3«2(T), 

w 4 (T) -> 6wa(T). (13) 

No matter what the form of gi(t), the final decay of the 
moments goes as v n (T) ~ 1/T. This is a consequence of 
the heavy weighting of <?i(i) near t = 0. 

C. Examples and Precedents 

Here we consider the intensity moments predicted by 
the above formalism for several forms of g\ (r) of experi- 
mental interest. Normalized variance predictions for five 
special cases are collected in Table [I] The first three 
of these cases are plotted vs exposure time in Fig. 0] 



These are larger for shorter exposures, and vanish for 
long exposures in the limit that the speckle is no longer 
visible. While dimensionless, these moments are not nor- 
malized in the sense that their values depend on the num- 
ber of speckles per pixel through (3. Hence we use a lower- 
case "«„(T)" for reduced moments defined in Eq. to 
contrast with the upper-case ll V 2 (T)" for the normalized 
second moment defined in Eq. (JSJ). Now the task is to 
compute an ensemble average over pixels i for quantities 
of form 



(10) 



(11) 
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FIG. 4: Comparison of [<?i(t)] 2 and Vz(T), for three types of 
scattering site dynamics, as labelled. According to prediction, 
Eq. JHj, the latter is a weighted average of the former. 



Note that the decay of V 2 (T) vs T is slower than that 
of gi(r) vs r; furthermore, the long-time decay is always 
V 2 {T) ~ 1/T, as explained by Eq. (JT3JI. This feature al- 
lows the characteristic time scale for the decay of <7i(r), 
which is equivalent to the characteristic broadening of 
the power spectrum, to be extracted even when it is a 
decade or more faster than the bandwidth of the camera. 

The only case for which we have analytically computed 
both second and third moments of the pixel signal distri- 
bution is for a Lorentzian spectrum or, equivalently, for 
an exponential field autocorrelation g\{r) = exp(— Tt). 



Invoking third and fourth order Siegert relations fl3| , and assuming that g\ (t) is even, we arrive at 

f T rh 

v 2 (T) = 2(3 -{l- t /T)\ gi {t)f, 

«3(T) = 3v 2 (T) + 12(3 2 J J -^(l-t 2 /T)g 1 (t 1 )g 1 (t 2 )g 1 (t 2 ~t 1 ), 

v 4 (T) = 6v 2 (T) + 0((3 2 ) 

I 
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gi(x = TT) 



V 2 (x = IT) 



cxp(— k) 
exp(- v / a) 
cxp(— x 2 ) 
\/6x / sinh \/6x 
V 6x 2 / sinh V 6x 2 



[e- 21 - 1 + 2x]/(2x 2 ) 
[(3 + + 4x)e- 2 ^ - 3 + 2x]/(2x 2 ) 
[e" 2 * 2 - 1 + v^zerf^x)]/^ 2 ) 
C(3)/x « 1.202/x, x > 1 
V / tt 4 /54/x « 1.343/x, x > 1 



TABLE I: Normalized variance predictions, computed from 
Eq. for various forms of the normalized electric field au- 
tocorrelation function, g\(j). The first example, <?i(t) = 
exp(— rr), corresponds to single-scattering from sites with 
diffusive dynamics and to DWS in backscattering for a sam- 
ple with ballistic dynamics. The second example, g\ (t) = 
exp(— yTV), corresponds to DWS in backscattering for a sam- 
ple with diffusive dynamics. The third example, gi(r) = 
exp[— (rr) 2 ], corresponds to single-scattering for a sample 
with ballistic dynamics. These three cases are plotted vs ex- 
posure time in Fig. 2] The fourth and fifth examples cor- 
respond to DWS in transmission for samples with diffuse 
and ballistic dynamics, respectively. These forms are too in- 
tractable to compute for general exposure times. 



This corresponds to single-scattering from a sample with 
diffusive dynamics and to DWS in backscattering from a 
sample with random ballistic dynamics. For the former, 
the linewidth or decay rate is T = Dq 2 where D is the 
diffusion coefficient and q is the magnitude of the scat- 
tering vector; for the latter, the decay rate isT« Att5v/X 
where Sv is the root-mean squared average random speed 
and A is the wavelength of light in the medium. For this 
example, the reduced second and third moments are 



-2x 



1 + 2x 



v 3 (T) = 3, 2 (T) + 6/3 2 (1 + a:)e ^3 1 + \ (14) 

where x — TT is the product of decay rate and exposure 
time, as per the notation in Table [I] 

The normalized variance for the special case of a 
Lorentzian spectrum, given in Eq. Ijl4ll . appeared nearly 
thirty years ago as Eq. (50) of Ref. 36]. It was subse- 
quently tested experimentally in Ref. |37| . This supports 
our theory of SVS, which seems both simpler and more 
general than that of Ref. j^| • Our approach applies for 
any form of 3i(t), not just for a Lorentzian spectrum, 
and it also accounts for any number of speckles per pixel. 
To our knowledge, Eqs. (|llb -c) have not previously ap- 
peared in the literature. 

The special case of a Lorentzian spectrum was also 
considered in Ref. |2|j , which is widely cited as a founding 
paper in the field of laser-speckle flowmetry. There the 
visibility of a speckle pattern is quantified by "speckle 
contrast" : 



K(T) = a T /(I), 



(15) 



quantity equals the square root of our reduced variance, 
K{T) = \A>2(r). The quoted result for a Lorentzian 
spectrum, Eq. (9) of Ref. [2jj and Eq. (13) of a more 
recent review |7J, would give V2(T) as the following un- 
weighted average of [31 (t)] 2 over the exposure interval: 



[ gi (t)] 2 dt/T = 



2x 



(16) 



where <jt is the standard deviation of the set of intensi- 
ties as measured over an exposure of duration T. This 



where x — TT as before. This conflicts with Eq. i|14|) 
here, and with Eq. (50) of Ref. 36], due to absence of 
the factors (3 and 2(1 - t/T) in Eq. ©. The latter mis- 
take of Ref. [2!| is that the variance is taken as a single 
integral of [31 (t)] 2 over the exposure window < r < T, 
rather than as a double integral where r = t\ — t% ranges 
over possible time differences within the window. The 
former mistake is that the value of (3, in effect, is taken 
as one; this is correct only if both the pixel size is in- 
finitesimal compared to speckle size and if just one po- 
larization mode is detected. A sampling of papers that 
invoke Eq. (9) of Ref. [2|j simultaneously match pixel size 
to speckle size but neglect an unknown visibility reduc- 
tion that results. The combined error introduced by the 
incorrect weighting and the neglect of j3 depend on de- 
tails of the experiment, but can easily exceed a factor of 
ten. Hence these issues may well contribute to the inabil- 
ity in the field of laser-speckle flowmetry to make repro- 
ducible quantitative connection between speckle visibility 
and blood flow speed. 



IV. DEMONSTRATION OF SVS 

In this section we both demonstrate the SVS technique 
and compare the experimental results with theoretical 
predictions of the previous section. 



A. Colloidal Particles 

Our first sample is the same opaque colloidal suspen- 
sion, probed by diffusely backscattered light with the 
same optical setup, as previously in Figs. 12131 For this 
sample, the speckles fluctuate due to diffusion of the par- 
ticles. Now we measure the second, third, and fourth 
moments of the distribution of pixel signals, Eq. and 
reduce the results to dimensionless form as per Eq. JSJ). 
This is done for many different exposure times T, with 
results shown vs T in Fig. [S] These reduced moments 
appear to approach a constant for short exposures, and 
to decay according to Eq. i|13fl as 1/T for long exposures. 

To compare with expectation, we first note that V2 (T) 
is a weighted average of /3[<?i(t)] 2 over the exposure in- 
terval, Eq. (|11|) . Therefore, we also include data for 
the latter as presented previously in Fig. [21 Recall 
that the functional form for the field autocorrelation 
is 3i (t) = exp(-VlV) with T = 6~/ 2 /t = 1623 s -1 . 
Evidently V2(T) and /3[3i(r)] 2 extrapolate to the same 
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FIG. 5: Reduced moments of the speckle pattern vs the ex- 
posure duration T, for light diffusely backscattering from an 
opaque colloidal suspension. Intensity autocorrelation data 
and fit from Fig. [5] are included for comparison. The solid 
curves for n = {2, 3} are generated numerically from the <?i(t) 
data in Fig. according to our theory of SVS, Eq. Hill . The 
dotted lines represent the expected long-exposure behavior, 
v n (T) ~ 1/T of Eq. The dashed curve is Eq. (JT^L the 

expectation for V2 (T) based on the formalism of Ref. 29) . 



FIG. 6: Reduced moments of the speckle pattern vs the ex- 
posure duration, T, for light diffusely transmitted through an 
opaque aqueous foam of thickness L = 1 cm. The solid curves 
represent numerical integration of the known field autocorre- 
lation function, gi(r) — y6Tr/ sinh \/6Tr, according to the 
prescription of Eq. li lt . The best fit is attained for /3 — 0.19 
and T = 121 s" 1 . 



value at short times, = 0.34. But while gi(r) decays 
more rapidly, ^(T) decays more slowly in accord with 
the heavy short-time weighting in the average across the 
exposure interval, Eq. I|ll|l . This qualitative agreement 
with expectation also can be made quantitative. Indeed, 
the two solid curves in Fig. are generated by numerical 
integration of <?i(t) data according to our SVS predic- 
tions of Eq. I|ll|) . As a check, the numerical prediction 
for V2(T) matches the analytic prediction given in Ta- 
ble The predictions for both V2(T) and vs(T) match 
the reduced moment data very well, with no adjustable 
parameters. 

Finally we compare with exp ectation based on the 
mistaken formalism of Ref. 29]. Introducing the cor- 
rect factor of (3 and taking the field autocorrelation as 
<?i(t) = exp(— VTV), the prediction for V2(T) would be 

/%[P,( i )Pd 1 /r = f 3 1 - (1 + 2 f S[1 "- 2 ^ ) , (it) 

where x = TT. This is plotted as a dashed curve for 
r = 1623 s _1 and (3 — 0.34, as known from the inten- 
sity autocorrelation data. Evidently, the formalism of 
Ref. (2i| does not correctly predict speckle variance. 

B. Foam 

As another example, we collect SVS data for light dif- 
fusely transmitted through an aqueous foam of thickness 
L = 1 cm [Gillette Foamy Regular] . The optical set-up is 
similar to the colloid experiments, except that the camera 
is moved opposite to the side upon which the laser light 
is incident. For this experiment, the speckles fluctuate 



due to sudden avalanche-like rearra ngem ents of bubbles 
within small localized sub- volumes |43|. Such dynam- 
ics are driven by the coarsening process, whereby small 
bubbles shrink and large bubbles grow in order to lower 
the total interfacial surface energy El]. Since the sam- 
ple is far from equilibrium, and the dynamics evolve with 
time, we restrict data collection to a narrow time window 
centered at 100 minutes after production. Here the aver- 
age bubble diameter is D ps 60 /im, the transport-mean 
free path is f « 3. 51?, and the average time between 
rearrangements at each scattering site is r Q ps 20 s [44^ . 
The volume of foam sampled by the collected photons 
is sufficiently great that the speckle pattern is in con- 
tinuous motion. The field correlation function takes the 
same form as for DWS in transmission from a sample 
of diffusing particles, gi{r) = %/ 6Tt / sinh V 6Tt. The 
first cumulant or initial decay rate is expected to be 
r = (L/1*) 2 /t ps (1 cm/0.02 cm) 2 /(20 s) = 125 s" 1 . 
This understanding is supported by both single-detector 
El El El EH and multi-speckle El El dynamic light 
scattering experiments. 

Our results for the second and third reduced moments 
of the distribution of pixel signals are shown in Fig. 
For short exposures both v^iT) and vs(T) approach a 
constant, from which we extrapolate to zero to find (3 = 
0.19. For longer exposures, the moments become smaller 
as the speckle pattern fluctuates more extensively during 
the exposure. To model this, we numerically integrate 
the field correlation function gi(r) = V 6Tt / sinh V 6Tt, 
according to the SVS prescription of Eqs. lfTT|) . Taking 
the first cumulant as T = 121 s _1 , close to expectation, 
we obtain a satisfactory fit to both v% (T) and v% (T) data 
as shown. 
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V. EXPERIMENTAL CONSIDERATIONS 

This final section provides guidance on the optimal 
design of an SVS experiment. Many of the issues, and 
the recommendations, are identical for other types of dy- 
namic light scattering experiment. Throughout we shall 
assume that statistics are not limited by lack of photons. 
In this case it is advantageous to double the laser power 
and to place a polarizer in front of the detector. While 
this does not change the average detected intensity, it 
does improve the contrast in intensity levels at the plane 
of the detector, and hence the signal-to-noise, since each 
polarization mode forms an independent speckle pattern. 
Therefore, throughout, we shall assume that polarized 
detection is employed. 



A. Optics 

First we consider the geometry of illumination and de- 
tection. Let a be the size of the region from which emerg- 
ing light is collected. For single-scattering experiments 
this could be controlled by the diameter of the incident 
beam or the length it travels within the sample. For 
multiple-scattering experiments it could be controlled by 
the beam diameter or the sample thickness. The value 
of a can also be affected by use of lenses or apertures 
between the sample and the detector. This is an impor- 
tant parameter because the angular size of the speckle, in 
the far field, is approximately A/a just as in a diffraction 
experiment. Thus, if the detector is located a distance 
d away from the source of the collected light, then the 
speckle size or spatial correlation length is approximately 
s = dX/a. 

Imagining that A/a is fixed, and that the light inten- 
sity can be adjusted at will, we now seek to optimize 
the distance d at which to place the detector. If d is 
too small, then the number N s of speckles at each pixel 
will be large and the intercept, or maximum contrast, 
(i = v 2 {0) = (I 2 )/(I) 2 - 1, will be small. The best case 
in terms of contrast is f3 — *■ 1 (or (3 — > 1/2 for unpolar- 
ized detection). In the opposite extreme, if d is too large 
then each speckle will span many pixels and the statistics 
of ensemble averaging will be poor. Overall, the figure 
of merit to be maximized is thus j3\fN s , the product of 
maximum contrast times the spread in number of speck- 
les per pixel. 

To find the optimal detector location by maximizing 
the figure of merit requires knowing [3 as a function of 
N s . We do this by Monte-Carlo simulation, calculat- 
ing the second intensity moment across a specified area 
for speckle patterns generated at random with the cor- 
rect statistical properties. Results for f3, as well as for 
the figure of merit (i\/N~ s , are plotted vs N s in Fig. [7| 
As expected, (3 — *■ 1 for small N s and /3 — * 1/N S for 
large N s . The figure of merit achieves a maximum where 
the speckle size nearly matches the pixel size, N s = 1. 
As Fig. \J\ demonstrates, an experiment is within about 
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FIG. 7: The intercept or maximum speckle contrast, /3 = 
v 2 (0) = {I 2 } /{I} 2 - 1, and the figure of merit fJ^/JH to be 
maximized in design of experiment, as a function of the num- 
ber N s of speckles per pixel. An experiment with polarized 
detection is within about ten percent of optimal if the inter- 
cept lies in the range 0.3 < (3 < 0.7. 



10% of optimum if the intercept lies within the range 
0.3 < (3 < 0.7 (or 0.15 < /3 < 0.35 for unpolarized de- 
tection). Thus a good strategy is to adjust the detector 
location until this criterion is met, keeping the illumina- 
tion optics fixed. It is well-known that pixel and speckle 
sizes should be matched, but to our knowledge specific 
guidelines in terms of the measurable intercept have not 
been published. 



B. Light Intensity 

Now we consider the optimal average intensity level, 
as controlled by choice of laser power. To beat photon- 
counting number fluctuations and dark-count subtraction 
error, this power should be as great as possible. How- 
ever, high power can result in clipping of signal for bright 
speckles that exceed the maximum grayscale level of the 
detector. This effect introduces an error whereby the 
measured intensity moments are shifted systematically 
to lower values. In the opposite extreme, for low laser 
power, the detected intensity levels are binned coarsely 
over too few grayscale levels. This effect introduces an 
error whereby the moments are shifted systematically to 
higher values. Two other effects can introduce systematic 
error at low laser power. One source is dark counts. For 
example, the pixels of our CCD camera report fluctuat- 
ing grayscale values of either 3 or 4, with a time average 
of 3.5, when there is no illumination. The other source 
of error is that the grayscale levels are reported at the 
lower edge of the bin, i.e. — 255 for an 8-bit camera 
like ours. For example, an actual signal level lying in the 
range 5 < S < 6 is reported as a grayscale level of 5. 

To investigate these effects we again turn to Monte- 
Carlo simulation. At first we restrict attention to an 
8-bit detector with a pixel size of three speckles. This 
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gives (3 w 0.3, and hence corresponds well with our col- 
loid experiments. In Fig. |S| we display results for the 
systematic error in the first four moments as a function 
of average intensity level. The top plot is for an ideal de- 
tector, with zero dark counts, with signal levels taken at 
the lower edge of the bins, {0, 1,2,..., 255}. Higher in- 
tensity levels are "clipped" to a value of 255. The fraction 
of pixels that must be clipped is plotted on the right-hand 
axis. At higher average intensity levels, where clipping 
occurs, the intensity moments fall below their correct val- 
ues. At lower average intensity levels, where digitization 
issues occur, the intensity moments rise above their cor- 
rect values. The middle plot shows that the latter can 
be mitigated to large extent by taking the signal level 
at the center of the bin. In other words, intensity mo- 
ments are much more accurate if an offset of 1/2 is added 
to each reported signal, so that possible levels are now 
{0.5, 1.5, 2.5, . . . , 255.5}. The bottom plot shows the ef- 
fect of dark counts, as simulated by randomly adding 3 
or 4 to the analog signal. This choice mimics the condi- 
tions of our colloid experiments. To mitigate both dark 
counts and lower-edge binning effect, we now subtract 
3.5 from each pixel value. Effectively, this introduces a 
±1 statistical error in pixel values, which broadens the 
distribution and causes higher than expected moments, 
as seen by comparison of Fig. |Hb-c. Under the operating 
conditions of our colloid experiment, (J) = 40, we use 
Fig. to estimate that the systematic error in our SVS 
data due to the combined effects of clipping, digitization, 
and dark-count effects is less than 0.25%. 

We now repeat the simulations for different numbers 
of speckles per pixel, assuming an 8-bit camera with zero 
dark counts. Plots of error vs average grayscale level are 
used to identify a safe operating range where the error in 
the variance is less than 0.1%. Recommended grayscale 
levels are shown as a function of intercept, (3, in Fig. 
Once an optimal value of (3 is achieved, for example by 
adjusting the detector location per the previous subsec- 
tion, the laser power should be adjusted according to this 
plot. Photon-counting and dark-count errors can be min- 
imized by operating at the upper end of the safe range. 



C. Normalization factor, j3 

The perceived contrast of the intensity levels in a 
speckle pattern is reduced progressively as the pixel size 
increases relative to speckle size. To eliminate this ef- 
fect, so that the remaining speckle contrast can serve 
as a quantitative probe of scattering site motion dur- 
ing the exposure window, the intercept (3 = t>2(0) = 
lim7^o((-f 2 )T/(-0 2 — 1) must be accurately determined. 
One approach, employed in our colloid and foam experi- 
ments above, is to collect data for many exposure times 
and to extrapolate the variance results to T = 0. This 
is satisfactory only if the dynamics are both stationary 
and sufficiently slow compared to the fastest speed of the 
camera. Obviously another approach is needed for sys- 
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FIG. 8: Simulated accuracy of SVS signals, due to errors in- 
troduced by an 8-bit digital camera, as a function of average 
grayscale level. For high average intensities, the brightest pix- 
els are clipped to a grayscale level of 255; the fraction of pixels 
that must be clipped is shown on the right-hand axis. The 
top plot shows results for moments computed directly from 
the returned grayscale levels {0, 1, 2, . . . , 255}. The middle 
shows how the accuracy is improved dramatically if pixel val- 
ues are offset by +1/2, which corresponds to the center of the 
grayscale bin. The bottom plot differs from the middle plot 
by inclusion of random dark counts, which are accounted for 
by subtraction. Note that the error scales in (b-c) are ten 
times smaller than the scale in (a). 



terns with fluctuating dynamics, where each individual 
exposure is to be analyzed in terms of scattering site mo- 
tion at that particular moment in time. This was the 
case for our first reported application of SVS, where we 
probed grain motion as a function of phase in a vibratory 
oscillation cycle |28| . There we had the luxury of being 
able to turn off the shaking and to measure the contrast 
of the static speckle pattern under absolutely identical 
illumination and detection conditions. 

Here we introduce an alternative method, whereby the 
value of can be eliminated from consideration alto- 
gether. The idea is to analyze not just one exposure, 
but rather some number m of successive exposures all 
of duration T. The first step is to find the variance for 
each of the exposures, and to average the results together, 
giving v 2 (T). The second step is to add together the m 
exposures pixel-by-pixel, and to compute the variance for 
the resulting "synthetic exposure" of duration mT, giv- 
ing V2 (mT) . These two variances depend on the value of 
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FIG. 9: Recommended average intensity levels for an 8-bit 
digital camera, as a function of the intercept — ^2(0); the 
corresponding pixel area is labelled in units of average speckle 
area. If the average intensity is too low, then error occurs due 
to digitizing the intensity crudely into a small number of bins. 
If the average intensity is too high, then error occurs due to 
clipping of the high-intensity tail of the distribution. The 
vertical bars indicate the range over which these effects cause 
less than a 0.1% error in the variance. For n-bit cameras, the 
lower limit does not change but the upper limit scales with 
the number of grayscale levels. 



[3, but their ratio does not: 

v 2 (mT) S™ T [l-t/{mT)][ gi {t)¥&t/m 



V2(T) 



/o T (l 



t/T)[ 9l {t)]Mt 



(18) 



The left-hand size is thus measured, and contact with 
scattering site motion is made by calculation of the right- 
hand side for the field autocorrelation of interest. The 
predicted forms in Table [I] can be used directly. For 
short exposures or slow dynamics, the variance ratios in 
Eq. i|18|) approach one. For long exposures or fast dy- 
namics, the variance ratios in Eq. (|18|) approach 1/m. 
As a specific example, the variance ratio for the case of a 
Lorentzian spectrum, g\{r) = cxp(— Tr), takes the form 



w 2 (mT) 

V2{T) 



— 1 + 2mx 



1 



1 + 2x)m 2 

12+2m+2ro 2 
15(l+m) X 



(19) 



-2m+3m" 2 



15(l+m) 



1 



2+2to+12to 2 
15(l+m) 



3-2m+3m 2 
15(l+ro) 



;,(20) 



where x = TT as before. The second line is a rational 
approximation that is correct to 0(x 3 ) and approaches 
1/m for long exposures. It can be inverted by solution 
of a quadratic equation. An additional advantage to this 
synthetic exposure method is that drift in laser power or 
detector gain, and CCD or CMOS noise that is correlated 
over successive exposures, are all automatically cancelled. 

This synthetic exposure variance ratio method is now 
illustrated for the same coarsening foam as in Fig. 
Here the foam is six hours old, and we employ an optical 
geometry whereby photons are both introduced and col- 
lected through the same 1 mm diameter aperture. This 
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FIG. 10: Pixel grayscale levels, variance ratios, and 
Lorentzian linewidths as a function of age for a coarsening 
foam. The offset along the time axis is to = 6 hrs, the lapse 
between sample production and the beginning of data collec- 
tion. Here the exposure duration is T = 10 ms. 



reduces the volume of foam sampled by the detected pho- 
tons, and ensures that only one rearrangement event is 
probed at a given time. Traditional DLS methods do not 
apply in this regime. An example rearrangement event 
is captured by SVS in Fig. ^| While the bubbles re- 
main in a fixed location, the speckle is nearly static and 
the variance ratios are nearly one. While the bubbles 
move, the speckle fluctuates and the variance ratios drop 
below one. According to the theory of DWS, the spec- 
trum is Lorentzian with linewidth T « AnSv/X, where 
Sv is the root-mean squared ballistic speed of the rear- 
ranging bubbles. Analyzing the variance ratio data using 
Eq. I|2UI) gives nearly identical linewidths for four differ- 
ent synthetic exposures, m = {2, 4, 8, 16}, as shown in the 
bottom plot of Fig.^H This good agreement further sup- 
ports our theoretical and experimental methods of SVS. 
It also demonstrates how rapidly varying dynamics may 
now be measured. 



Acknowledgments 

We thank M. Giglio, P.-A. Lcmieux, R.P. Ojha, P.N. 
Pusey, and T. Usher for helpful discussions. This mate- 
rial is based upon work supported by NSF under grants 
DMR-0305106 and PHY-0320752, and by NASA under 
Microgravity Fluid Physics grant NAG3-2481. 



12 



[1] H. Z. Cummins and E. R. Pike, eds., Photon Correla- 
tion and Light-Beating Spectroscopy (Plenum Press, NY, 
1974). 

[2] B. J. Berne and R. Pecora, Dynamic Light Scattering: 

With Applications to Chemistry, Biology, and Physics 

(Dover Publications, NY, 2000). 
[3] B. Chu, Laser Light Scattering, Basic Principles and 

Practice (Academic Press, NY, 1991). 
[4] W. Brown, ed., Dynamic Light Scattering: The Method 

and Some Applications (Claredon Press, Oxford, 1993). 
[5] A. P. Shepherd and P. A. Oberg, Laser-Doppler Blood 

Flowmetry, vol. 107 of Developments in Cardiovascular 

Medicine (Kluwer Academic Publishers, Boston, 1990). 
[6] Y. Aizu and T. Asakura, Opt. Laser Tech. 23, 205 (1991). 
[7] J. D. Briers, Physiol. Meas. 22, R35 (2001). 
[8] G. Maret and P. E. Wolf, Z. Phys. B 65, 409 (1987). 
[9] D. J. Pine, D. A. Weitz, P. M. Chaikin, and E. Her- 

bolzheimer, Phys. Rev. Lett. 60, 1134 (1988). 
[10] D. A. Weitz and D. J. Pine, in Dynamic Light Scattering, 

edited by W. Brown (Claredon, Oxford, 1993), p. 652. 
[11] G. Maret, Curr. Op. Coll. Int. Sci. 2, 251 (1997). 
[12] J. D. Briers, J. Opt. Soc. Am. A 13, 345 (1996). 
[13] P.-A. Lemieux and D. J. Durian, J. Opt. Soc. Am. A 16, 

1651 (1999). 

[14] W. C. K. Poon, Curr. Op. Coll. Int. Sci. 3, 593 (1998). 
[15] L. Cipelletti and L. Ramos, J. Phys.-Cond. Matt. 17, 
R253 (2005). 

[16] H. M. Jaeger, S. R. Nagel, and R. P. Behringer, Rev. 

Mod. Phys. 68, 1259 (1996). 
[17] A. Wong and P. Wiltzius, Rev. Sci. Inst. 64, 2547 (1993). 
[18] S. Kirsch, V. Frenz, W. Schartl, E. Bartsch, and 

H. Sillescu, J. Chem. Phys. 104, 1758 (1996). 
[19] L. Cipelletti and D. A. Weitz, Rev. Sci. Inst. 70, 3214 

(1999). 

[20] D. Lumma, L. B. Lurio, S. G. J. Mochrie, and M. Sutton, 

Rev. Sci. Inst. 71, 3274 (2000). 
[21] V. Viasnoff, F. Lequeux, and D. J. Pine, Rev. Sci. Inst. 

73, 2336 (2002). 
[22] J. Xu, X. Dong, L. F. Zhang, Y. G. Jiang, and L. W. 

Zhou, Rev. Sci. Inst. 73, 3575 (2002). 
[23] L. Cipelletti, H. Bissig, V. Trappe, P. Ballesta, and 

S. Mazoyer, J. Phys.-Cond. Matt. 15, S257 (2003). 
[24] T. Seydel, A. Madsen, M. Sprung, M. Tolan, G. Grubel, 

and W. Press, Rev. Sci. Inst. 74, 4033 (2003). 
[25] K. N. Pham, S. U. Egelhaaf, A. Moussaid, and P. N. 

Pusey, Rev. Sci. Inst. 75, 2419 (2004). 
[26] P. Falus, M. A. Borthwick, and S. G. J. Mochrie, Rev. 

Sci. Inst. 75, 4383 (2004). 



[27] P. K. Dixon and D. J. Durian, Bull. Am. Phys. Soc. 46, 
932 (2001). 

[28] P. K. Dixon and D. J. Durian, Phys. Rev. Lett. 90, 
184302 (2003). 

[29] A. F. Fercher and J. D. Briers, Opt. Comm. 37, 326 
(1981). 

[30] H. C. van de Hulst, Light scattering by small particles 

(Dover, New York, 1957). 
[31] A. A. Middleton and D. S. Fisher, Phys. Rev. B 43, 5934 

(1991). 

[32] D. J. Durian, Phys. Rev. E 51, 3350 (1995). 

[33] F. C. MacKintosh, J. X. Zhu, D. J. Pine, and D. A. Weitz, 

Phys. Rev. B 40, 9342 (1989). 
[34] M. U. Vera, P.-A. Lemieux, and D. J. Durian, J. Opt. 

Soc. Am. A 14, 2800 (1997). 
[35] P.-A. Lemieux, M. U. Vera, and D. J. Durian, Phys. Rev. 

E 57, 4498 (1998). 
[36] E. Jakeman and E. R. Pike, J. Phys. A (Proc. Phys. Soc.) 

1, 128 (1968). 

[37] E. Jakeman, C. J. Olivier, and E. R. Pike, J. Phys. A 

(Proc. Phys. Soc.) 1, 406 (1968). 
[38] A. K. Dunn, T. Bolay, M. A. Moskowitz, and D. A. Boas, 

J. Cer. Blood Flow Metab. 21, 195 (2001). 
[39] T. Durduran, M. G. Burnett, G. Q. Yu, C. Zhou, D. Fu- 

ruya, A. G. Yodh, J. A. Detre, and J. H. Greenberg, J. 

Cer. Blood Flow Metab. 24, 518 (2004). 
[40] B. Weber, C. Burger, M. T. Wyss, G. K. von Schulthess, 

F. Scheffold, and A. Buck, Eur. J. Neurosci. 20, 2664 

(2004). 

[41] S. Yuan, A. Devor, D. A. Boas, and A. K. Dunn, App. 
Opt. 44, 1823 (2005). 

[42] J. W. Goodman, in Laser Speckle and Related Phenom- 
ena, edited by J. Dainty (Springer- Verlag, Berlin, 1984), 
vol. 9 of Topics in Applied Physics, 2nd ed. 

[43] D. J. Durian, D. A. Weitz, and D. J. Pine, Science 252, 
686 (1991). 

[44] D. J. Durian, D. A. Weitz, and D. J. Pine, Phys. Rev. A 

44, R7902 (1991). 
[45] J. C. Earnshaw and A. H. Jaafar, Phys. Rev. E 49, 5408 

(1994). 

[46] A. D. Copal and D. J. Durian, J. Coll. I. Sci. 213, 169 
(1999). 

[47] S. Cohen-Addad and R. Hohler, Phys. Rev. Lett. 86, 
4700 (2001). 

[48] P. Mayer, H. Bissig, L. Berthier, L. Cipelletti, J. P. Gar- 
rahan, P. Sollich, and V. Trappe, Phys. Rev. Lett. 93, 
115701 (2004). 



